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The analysis of Balmcr-dominated optical spectra from non-radiative (adiabatic) SNRs has shown 
that the ratio of the electron to proton temperature at the blast wave is close to unity at vs ^400 
kms -1 , but declines sharply down to the minimum value of m e /m p dictated by the jump conditions 
at shock speeds exceeding 2000 kms -1 . We propose a physical model for the heating of electrons and 
ions in non-cosmic ray dominated, strong shocks (vs > 400 km s -1 ) wherein the electrons are heated 
by lower hybrid waves immediately ahead of the shock front. These waves arise naturally from the 
cosmic ray pressure gradient upstream from the shock. Our model predicts a nearly constant level of 
electron heating over a wide range of shock speeds, producing a relationship (T e /T p )o oc Vg 2 (oc M~ 2 ) 
that is fully consistent with the observations. 

Subject headings: ISM: kinematics and dynamics, shock waves, plasmas, ISM: cosmic rays, supernova 
remnants 



1. INTRODUCTION 

The discovery of collisionless shock waves in the solar 
wind in the 1960s ushered in a new era in the physics of 
space plasmas (see Tidman & Krall 1971, Sagdeev 1979 
and Kennel 1985 for reviews and references). Due to the 
low density (n<l cm -3 ) of the interplanetary medium, 
the jump in hydrodynamical quantities is produced not 
by Coulomb collisions but by collective plasma processes 
such as electromagnetic waves and turbulence. However, 
despite the availability of extensive in situ observations 
of solar wind shocks and the expenditure of consider- 
able effort in theoretical modelling of these structures, 
a detailed understanding of the processes at the shock 
transition responsible for partitioning the shock energy 
between different charged particle species has been slow 
to emerge. The problem is far more acute for inter- 
stellar shocks, where in situ measurements of the the 
shock structure are unavailable and the very high Mach 
numbers (^30—200) make numerical simulations of these 
structures extremely difficult or impossible. 

The optical emission generated by non-radiative super- 
nova remnants (SNRs) in partially neutral gas provides 
a valuable diagnostic tool for probing the heating pro- 
cesses in collisionless shocks. The optical spectra of these 
SNRs (which lose a negligible fraction of their energy to 
radiation) are dominated by Balmer line emission, pro- 
duced by collisional excitation when neutral hydrogen is 
overrun by the blast wave (Chevalier & Raymond 1978, 
Bychkov & Lebedev 1979). Each emission line consists 
of two components: (1) a narrow velocity component 
produced when cold, ambient H I overrun by the shock 
is excited by electron and proton collisions, and (2) a 
broad velocity component produced when fast neutrals 
created by postshock charge exchange are collisionally 
excited (Chevalier, Kirshner & Raymond 1980). An ex- 
ample of a Balmer-dominated shock spectrum from the 
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Galactic SNR RCW 86 is shown in Fig. 1. The opti- 
cal emission arises in a very thin (< 10 16 cm) ionization 
zone, thin enough so that the protons transformed into 
hot neutrals have had little time to equilibrate with elec- 
trons and other ions. This makes the measured width of 
the broad Balmer line directly proportional to the pro- 
ton temperature set by collisionless heating at the shock 
front. The broad-to-narrow flux ratio, Ib/In on the 
other hand, is sensitive to both the degree of electron- 
proton temperature equilibration at the shock front (i.e., 
(T e /T p )o ) and the shock velocity, vs (Chevalier, Kirsh- 
ner & Raymond 1980, Smith et al. 1991). The ratio also 
depends (though less sensitively) on the preshock neu- 
tral fraction. The broad component width and 1b /In 
of an observed Balmer-dominated shock can be modeled 
with numerical shock codes to simultaneously estimate 
v s and (T e /T p ) (Smith et al. 1991, Ghavamian 1999, 
Ghavamian et al. 2001, Ghavamian et al. 2002). 

In this paper we draw together observed values of 
(T e /T p )oand v$ measured in five Balmer-dominated 
SNRs, including several previously unpublished measure- 
ments from the SNR RCW 86. We then propose a phys- 
ical model of electron heating by lower hybrid waves in 
a cosmic-ray precursor that obeys the observed relation- 
ship between (T e /T p )o and vs ■ We conclude by explor- 
ing some consequences of this interpretation, the con- 
nection to other observables and applicability to other 
collisionless shock situations. 

2. RELATIONSHIP BETWEEN EQUILIBRATION AND 
SHOCK SPEED 

In Fig. 2 we plot 11 measurements of (T e /T p )o and 
vs , obtained from longslit spectra of four Galactic rem- 
nants: the Cygnus Loop (one position from the north- 
east), RCW 86 (positions from the southwestern, north- 
ern and eastern limbs), Tycho's SNR ('Knot g' from the 
eastern rim), SN 1006 (one position from the northwest- 
ern rim) and one remnant in the Large Magellanic Cloud 
DEM L 71 (five positions from the full circumference 
of the shell). In the first four remnants, the shock pa- 
rameters have been estimated via longslit spectroscopy. 
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Fig. 1. — An example of the optical spectrum of a Balmcr- 
dominatcd shock, showing the broad and narrow Ha lines char- 
acteristic of non-radiative shocks in partially neutral gas. This 
spectrum, originally presented by Sollcrman et al. (2003), was 
obtained from the southwestern rim of the Galactic SNR RCW 86, 
with high enough spectral resolution (~10 kms ) to resolve the 
broad (~ 500 kms -1 FWHM) and narrow (~30 kms -1 FWHM) 
Ha lines. The night sky OH lines (indicated by the © symbol) 
have been left in to demonstrate their relatively narrower widths 
compared to the Ha lines. The broad Ha width and ratio of the 
broad to narrow Ha flux for these types of shocks were used to 
produce the relationship shown in Fig. 2. 

First, wc used the broad component Ha width to con- 
strain the range of shock speeds for each Balmer fila- 
ment in the limits of minimal ((T e /T p )o = m e /m p ) and 
full ((T e /T p )o = 1) equilibration. Next, we fine tuned 
(T e /T p ) and v$ by using shock models to match the 
observed broad to narrow flux ratio (Ghavamian et al. 
2001; 2002). 

In the case of RCW 86 we have added two more 
data points, indicated by dashes in Fig. 2, derived from 
previously unpublished optical spectra of Ghavamian 
(1999). These were obtained from the eastern rim 
{vFWHM(Ra) = 640±35 kms" 1 ; I B /I N =1.0±0.2) and 
the northern rim (vpwHM(tlot) = 325±10 kms -1 ; 
I B /I N = 1.06±0.1) of RCW 86. In the fifth SNR of our 
sample, DEM L 71 (Ghavamian et al. 2003), we used 
the broad component Ha width only to constrain the 
range of shock speeds and the proton temperature T p , 
due to the anomalously low Ib/In values in this SNR 
(see below) . We then combined this information with T e 
measured from Chandra X-ray spectra of the blast wave 
to obtain (T e /T p ) (Rakowski, Ghavamian & Hughes 
2003; Rakowski 2005). 

Although the equilibrations were obtained from differ- 
ent SNRs in a range of environments and distances, they 
show a clear trend of decreasing (T e /T p )o with v$ ■ The 
plot in Fig. 2 suggests that a transitional shock speed 
exists around 400 kms -1 , below which collisionless pro- 
cesses promptly equilibrate the electrons and protons at 
the shock front. Above that speed (T e /T p )o rapidly de- 
clines, eventually reaching values consistent with mass- 
proportional heating. This decline appears to follow 



(T e /T p )o oc Vg 2 . We show that this behavior is pre- 
dicted by a physical model of electron heating, where 
the electrons are heated in the cosmic ray precursor to 
a level that is nearly independent of shock speed before 
acquiring the mass-proportional increment at the shock 
front. The protons, on the other hand, receive only mass- 
proportional heating, resulting in the inverse squared de- 
pendence of equilibration on shock speed. Note that 
while our proposed model relies on the acceleration of 
cosmic rays and the existence of a cosmic ray precursor 
to heat the electrons, it does not require the energetics 
of the shock to be dominated by the cosmic rays. 

In physical models the strength of the (assumed quasi- 
perpendicular) shock is characterized not by the shock 
speed, but rather the magnetosonic Mach number Ms 
(= vs/vms, where v M s = (c| + v^) 1 / 2 is the magne- 
tosonic speed, cs is the sound speed (= \J an d va 

(= B/y/An pi) the Alfven speed of the prcshock gas). 
The preshock temperature, ion density and magnetic 
field strength are not strongly constrained in the ob- 
served Balmer-dominated shocks. In particular, Ms is 
most sensitive to the choice of preshock magnetic field 
due to the dependence of the Alfven speed on B 2 . How- 
ever, assuming standard values for the warm neutral ISM 
- preshock temperature of 10,000 K, density of 1 cm~ 3 , 
magnetic field strength of 3 fxG and 50% preshock ion- 
ization, the magnetosonic speed is then approximately 
13 kms -1 . The corresponding values of Ms are marked 
at the top of Fig. 2. 

3. LOWER HYBRID WAVE HEATING MODEL 

The constant electron heating with shock velocity (giv- 
ing (T e /T p )o oc Wj ) suggests a process occurring in 
the preshock medium rather than the shock front it- 
self. Waves in the preshock medium can be excited 
by shock reflected ions which gyrate around the field 
lines before returning to the shock itself, or by a cos- 
mic ray precursor. Cargill & Papadopoulos (1988) per- 
formed hybrid simulations of the first possibility, wherein 
shock reflected ions generate Langmuir and ion-acoustic 
waves ahead of the shock, which then heat electrons as 
they are damped. They predicted a temperature ratio 
(T e /T p )o ^0.2 nearly independent of shock velocity, a re- 
sult which clearly disagrees with the observational data 
in Fig. 2. Laming (2001a, b) modeled the excitation of 
lower hybrid waves by shock reflected ions, following the 
suggestion of McClements et al. (1997) that lower-hybrid 
waves could stream away from the shock with a group ve- 
locity equal to the shock velocity. In this way, waves stay 
in contact with the shock for arbitrarily long periods of 
time and hence can grow to large intensities despite low 
intrinsic growth rates. The electron heating predicted 
by such a model (Vink & Laming 2003) depends on the 
product of the maximum electron energy and the fraction 
of electrons which are resonant with the waves that can 
be accelerated. This gives (Te/T^ocx^ 1 xcxp(-M 2 ), 
where M is the shock Mach number. 

Here we explore another approach. Cosmic rays have 
been long known to generate waves upstream of shocks; 
the generation of Alfven waves is an intrinsic part of 
cosmic ray acceleration models. Drury & Falle (1986) 
showed that the negative gradient of cosmic pressure 
with distance ahead of the shock can generate sound 
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waves via the "Drury instability" . This is the mecha- 
nism which smoothes out the hydrodynamical jump in 
cosmic ray modified shocks. Mikhailovskii (1992) gives a 
corresponding expression for the excitation of magneto- 
acoustic waves. We argue that the high frequency exten- 
sion of such processes by the cosmic ray pressure gradi- 
ent would be the excitation of lower- hybrid waves. The 
frequencies of these waves would lie between the gyrofre- 
quencies of the electrons and the protons, with electron 
heating occurring as the waves damp along magnetic field 
lines. Without calculating the growth rate of the lower 
hybrid waves explicitly, we can estimate the magnitude 
of the electron heating. The parallel diffusion coefficient 
for electrons in lower hybrid wave turbulence is (derived 
from Equations 10.83 and 10.93 of Melrose (1986); see 
also Begelman & Chiueh 1988) 



4 V m e J k 2 . 
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where u> = ^n e Q p is the lower hybrid wave frequency, 
the geometric mean of the electron and proton cyclotron 
frequencies (Q e _ p = eB /m e ^ p c), and fcii and k±_ are 
wavevectors parallel and perpendicular to the magnetic 
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field. Taking SE = B f^fj ' ^ (Karney 1978), 
kf\/k± = m e/ m p, an d the time spent by an electron in- 
side the cosmic ray precursor t ~ n/vg where k is the 
cosmic ray diffusion coefficient, we find 
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If the lower hybrid wave group velocity du>/dk± — vg, 
then the phase velocity uj/k± = 2vg when k^/k 2 ^ = 

m e /m p (Laming 2001a), and the electron heating de- 
pends on a few constants times the product £l e K. Thus 
for a Bohm-like cosmic ray diffusion coefficient n cx 1/B 
and independent of vg , we have an electron heating pro- 
cess independent of both the shock velocity and the up- 
stream magnetic field, and (T e /T p )o oc Vg as required. 
This estimate neglects the fact that only electrons with 
velocity greater than v e ~ w/fcii may be heated by the 
waves. However, for the likely depth of a cosmic ray 
precursor (estimated below), the electrons will have suf- 
ficient time to at least partially collisionally equilibrate 
amongst themselves before crossing the shock, allowing 
a much larger fraction of preshock electrons to interact 
with the waves. The time spent by an electron inside a 
reflected ion precursor, as opposed to a cosmic ray pre- 
cursor, would be t ~ d/vg ~ l/f2 p , where the precursor 
depth d is approximately the gyroradius of a cosmic ray 
proton. Substituting this relation into equation 2 gives 
(T e /T p )o independent of vs , similar to the behavior pre- 
dicted by Cargill & Papadopoulos (1988), although these 
authors consider different waves. In this case the shallow- 
ness of the reflected ion precursor prevents the collisional 
redistribution of energy amongst the electrons. 




1000 2000 3000 

Shock Velocity (km s _1 ) 

Fig. 2. — The electron to proton temperature ratio at the 
shock front as a function of shock velocity for 5 Balmer-dominated 
SNRs. Magnetosonic Mach numbers (Ms) appropriate for typical 
ISM conditions are indicated along the top axis. The data shown 
here were measured from Balmer-dominated shocks in the Cygnus 
Loop, RCW 86, Tycho's SNR (Ghavamian et al. 2001), SN 1006 
(Ghavamian et al. 2002) and DEM L71 (Rakowski et al. 2003, 
Rakowski 2005). The dashed error bars for RCW 86 mark previ- 
ously unpublished results. Below 400 km s — 1 (Ms ~ 30) the data 
are consistent with (T e /T p )o =1. The prediction of the proposed 
lower hybrid wave heating mechanism in the cosmic-ray precursor, 
(T e /T p )o ocVg 2 (oc Mg 2 ), is shown for v s > 400 krns" 1 . 

If a constant kinetic energy is imparted to electrons 
by lower hybrid waves, then for shock speeds of 400 
kms -1 and above (c.f. Fig. 2) the observations require 
/2 ~ 4 x 10~ 10 ergs (~ 0.3 keV) immediately be- 
hind the shock front. In that case, we infer £l e n ~ 7x 10 21 
cm 2 s~ 2 , and k ~ 4 x 10 14 /B cm 2 s _1 . For typical 
magnetic fields B ~3 [iG, k is several orders of magni- 
tude below that inferred for the undisturbed interstellar 
medium, but entirely consistent with values estimated 
for the solar wind or interplanetary medium. We em- 
phasize that this estimate is a lower limit on k. This 
estimate neglects the electron-electron collisional equili- 
bration necessary to bring more electrons into resonance 
with the turbulence. Furthermore, lower hybrid waves 
can only be excited in the cosmic ray precursor if the cos- 
mic ray pressure gradient is sufficiently strong. The re- 
quired gradient scale lengths are 1/L > 8cs /3k for sound 
waves (Drury & Falle 1986), and 1/L > VA/cr g ~ va/k 
for magneto-acoustic waves (Mikhailovskii 1992), where 
r g is the gyroradius of a cosmic ray proton. An up- 
per limit on k will come from the requirement that the 
neutral H survive against electron impact ionization in 
the precursor to reach the shock front. Numerically, 

k < 10 24 (w s /1000 km s" 1 ) /n e cm s . Further con- 
straints on the cosmic ray diffusion coefficient and pres- 
sure will be obtained from calculation of the lower hybrid 
wave growth rate required to sustain the electron heat- 
ing. This growth rate must be significantly larger than 
the charge exchange frequency of the partially neutral 
gas. Under typical ISM conditions ^Cl e D, p ~ 1 Hz and 



4 



una (cr cx v) ~ 10~ 8 Hz, so the condition for electron heat- 
ing is nearly always satisfied. 

The electron heating by lower-hybrid waves is inher- 
ently anisotropic. If some of this anisotropy survives 
the electron-electron collision equilibration, a polariza- 
tion signal may be present in the narrow component of 
Ha (see e.g. Laming 1990), which might provide more 
insight into the electron heating process. 

4. THE WIDTH OF THE Ha NARROW COMPONENT 

The same cosmic ray precursor should also generate 
lower frequency waves, which for a quasi-perpendicular 
shock will include magneto-acoustic waves. Below the 
ion-neutral charge exchange frequency these waves are 
not effectively damped (Drury et al. 1996), and may re- 
veal themselves via broadening of the narrow Ha compo- 
nent (Smith, Raymond, & Laming 1994). Those authors 
suggested that waves in the cosmic ray precursor actu- 
ally heat the preshock gas to 30,000 - 40,000 K and that 
the observed linewidth is thermal in nature. An upper 
limit to the cosmic ray diffusion coefficient of ~ 10 24 /n e 
cm 2 s _1 then results from the constraint that sufficient 
neutral H survive passage through the precursor to en- 
counter the shock. However no precise heating mecha- 
nism was specified, and we now suggest that the large 
width of the narrow component Ha line in the observed 
SNRs is not of thermal origin, but rather due to the 
motion of protons in the lowest frequency waves of the 
magneto-acoustic spectrum. These waves lie below the 
charge exchange frequency, allowing a coherent oscilla- 
tion of preshock neutrals and preshock protons. In this 
case, Sv ~ vaSB/B for Alfven waves (va should be re- 
placed by vms f° r magneto- acoustic waves). The Alfven 
speed upstream from the shock is in the range 1-10 km 
s _1 , so an observed broadening Sv/va > 1 also implies 
SB/B > 1. Such a magnetic field amplification has al- 
ready been proposed in the context of cosmic ray ac- 
celeration in shocks (Lucek & Bell 2000, Bell & Lucek 
2001). A successful model of precursor ions will need to 
reproduce the existing measurements of narrow compo- 
nent Ha line widths (Smith, Raymond & Laming 1994; 
Hester et al. 1994; Sollerman ct al. 2003) which indi- 
cate that they remain relatively constant (30-50 kms -1 ) 
over a wide range in shock speed (180-3000 kms -1 ). We 
defer this topic to future work. 

5. DISCUSSION AND SUMMARY 

The plot of (T e /T p )oin Fig. 2 is qualitatively similar 
to results of a survey of electron heating at solar wind 
shocks (Schwartz et al. 1988). The main difference is 
that T e <~ Ti is obtained for Alfvenic Mach numbers 



< 5 in the solar wind, whereas in Fig. 2 this occurs 
for magnetosonic Mach numbers around 20-30, for as- 
sumed preshock magnetic fields of 3^iG. A preshock mag- 
netic field amplification in SNRs, of a similar amount 
required to reproduce the width of the Ha narrow com- 
ponent (SB/B ~ 5-10), would shift the {T e /T p ) ratios 
in Fig. 2 onto Mach numbers similar to those in the solar 
wind. This may support our interpretation. Collisionlcss 
shocks play a dominant role in heating the intergalactic 
medium (IGM) during large scale structure formation 
(Ryu et al. 2003; Yoshida et al. 2005; Kang et al. 
2005). If the preshock magnetic field is amplified as de- 
scribed above, then the relationship (T e /T p ) cx Mg 2 
may also be applicable to modeling the ionization and 
temperature structure of the IGM. 

The maximum energy attainable by cosmic rays 
in Balmer-dominated shocks is limited by ion-neutral 
damping to -0.1-1 TeV (Drury et al. 1996). It is 
possible that the shocks with stronger cosmic ray ac- 
celeration generate sufficient electron heating in their 
precursors that no neutrals survive to cross the shock. 
The anomalously low Ig /In ratio observed in DEM L71 
(Ghavamian et al. 2003) may indicate that this is an 
intermediate case. Electron heating in a shock precur- 
sor will produce added narrow component Ha. However, 
since the broad component can only arise by charge ex- 
change with shocked protons, the neutral H must pene- 
trate either the shock or the reflected ion precursor. We 
emphasize that of the points plotted in Fig. 2 from Ha 
Ib/In measurements, the electron temperatures from 
RCW 86 (Vink et al. 2006), the Cygnus Loop (Lev- 
enson et al. 2002), SN 1006 (Laming et al. 1996, Vink 
et al. 2003) and Tycho (Hwang et al. 2002; Warren 
et al. 2005) have been independently corroborated by 
measurements in other wavebands. 

Several lingering questions remain from Fig. 2: what 
mechanism causes the prompt electron-proton equilibra- 
tion below 400 kms -1 , and why does the decline in 
(T e /T p ) begin at that shock speed? Even more im- 
portantly, does the inverse square relationship between 
equilibration and shock speed/Mach number also hold 
for collisionless shocks in fully ionized gas? These issues 
will be addressed in future work. 
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